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Abstract 

We calculate, in the standard Bogoliubov approximation, the ground state 
energy of the spinor BEC with hyperfine spin / = 1 where the two-body re- 
pulsive hard-core and spin exchange interactions are both included. The cou- 
pling constants characterized these two competing interactions are expressed 
in terms of the corresponding s-wave scattering lengths using second-order 
perturbation methods. We show that the ultraviolet divergence arising in the 

ground state energy corrections can be exactly eliminated. 
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Recently, Stamper-Kurn et al. |]] have successfully cooled 23 Na atoms using an optical 
dipole trap and achieved BEC. In their experiment, a multi-component BEC which is char- 
acterized by the three hyperfme spin states |/ = 1, m./ = ±1, 0) has been observed. This 
has opened an interesting possibility to explore the multi-component BEC with complicated 
internal spin dynamics in which not only the global U (1) symmetry but also the rotational 
SO (3) symmetry in spin space are involved P|§. 

An important feature of the spinor condensate is that, in addition to the repulsive bi- 
nary hard-core collisions which give rise to the density- density interaction, atoms in the 
condensates can also couple to each other via the spin exchange interaction. Assuming that 
the interaction for each spin exchange channel is again characterized by zero-range delta- 
potential scattering, one thus obtains the interacting term S • S where S is the spin density 
operator. The competition between these two interactions thus lead to an intriguing scenario 
of the spin dynamics of the spinor BEC which is characterized by a complex ground state 
structure @-§j. 

The question arise now that how these two-body interactions alter the dynamical prop- 
erties of the condensates. It is known that in a weakly interacting Bose condensed system, 
the two-body interactions play a crucial role in determine the low temperature properties of 
the systems, which will modify the ground state of the many-particle systems and cause a 
depletion of the condensate fraction even at the zero temperature. Moreover, an divergence 
could possibly appear when we calculate the ground state energy in the standard Bogoliubov 
approximation. This divergence is due to the naive assumption of a constant matrix ele- 
ment of binary interaction irrespective of the relative momenta of the interacting particles. 
A well illustrated example is the ultraviolet divergence occurring in the calculated ground 
state energy of the one-component BEC where the two-body interaction is described by 
the repulsive hard-core collisions with a momentum- independent coupling constant [[J. To 
eliminate such a divergence and gain more insights into the ground state properties of the 
condensates, one has to calculate the s-wave scattering length at least to the second order in 
coupling constant. By expanding the coupling constant in powers of the s-wave scattering 
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length the previously mentioned divergent ground state energy can be rendered finite. This 
is expected since in a physically sensible theory the ground state energy must assume a finite 
value when expressed in terms of physically measurable quantities. In other generalized Bose 
condensed systems such as the spinor BEC ultraviolet divergences of the same sort could 
also appear. It is therefore quite essential to verify if a similar procedure could completely 
removed these divergences, and in this paper we address this issue in details for the / = 1 
spinor BEC in the presence of a constant magnetic field. 

Consider an assembly of homogeneous dilute Bose gas with hyperfine spin / = 1. The 
natural basis set to characterize such a system is the hyperfine spin states |m/ = ±1,0} . 
However, in view of the special symmetrical forms of the S = 1 spin matrix representations 
, one may adopt the basis set {\x) , \y) , \z)} which is defined as the eigenstates of the a-th 
component of the spin operator with eigenvalue 0, i.e., S a \a) — (a = x, y, z) such that the 
matrix elements for the a-th spin component is given by (7 \S a \ (5) = ie al 3 y where e a/ g 7 is the 
Levi-Civita tensor. This representation enables us to relate the matrix elements of S = 1 
spin operators to those of the space rotation, allowing the order parameter to behave as a 
vector under spin space rotation. 

For the / = 1 spinor BEC, the bosonic atomic field can be described by the multi- 
component field operator with components ipa ( r ) i a — X ,V, Z ), an d thus the density of 
the particle number and spin can be written as n = ip^ipa, and S a = ipaS a ip/3 = —iea^iptip^ 
respectively. Note that we have used the summation convention over the indices of com- 
ponent a, (3, ■■■ throughout this paper. Now, without loss of generality, the Hamiltonian 
density can be constructed in the presence of a constant magnetic field B pointing to the 
^-direction: 

h = -^l^a + \g n n 2 + ^ 9s s • s - n • s (h = 1) (1) 

where Q =Qz = g M B (g M : gyromagnetic ratio) is the Larmor frequency in a vectorial 

notation. Expanding n and S a in terms of the field operators, Eq.(|T]) can be expressed 
as 
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where the two new coupling constants are given by g\ — g n + g s , g% — —g s . According to 
the recent spectroscopic experiment by Abraham et al. 0, it is conceivable in general that 
#2 is comparable to g\ in magnitude and can be either positive or negative. It is known 
that the positive g2 implies the ferromagnetic coupling while the negative one implies the 
antiferromagnetic coupling for the spin exchange interaction. 

Since the system is homogeneous, the field operator can be expanded in terms of creation 
and annihilation operators characterized by momentum k 

^(r) = 47E«^ kr , (3) 

vv k 

where V denotes the volume of the system. Accordingly, the Hamiltonian in momentum 
space now reads as 

H = Hq + -ffmag + H int (4) 



where 



Hq = 53 e kal,k a ",k (5) 
k 

(6) 



IT - 3l \ - ft 

tl mt - — 2^ a /3,k 4 a a,k 3 a a,k2 a /3,ki 

ki+k2=k3+k4 

_l_ 32 ST^ t t (n\ 

2V a /3,k 4 a /3,k 3 a Q;,k2 a a,ki {' ) 

ki+k2=k 3 +k4 

In the ground state, most particles occupy the k = states. As a result, the scatter- 
ing between two nonzero-momentum states can be ignored and the interacting part of the 
Hamiltonian can be replaced by 



mt ~ 2V 
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Before calculating the s-wave scattering lengths to the second order, a couple of re- 
marks are in orders. First of all, for the sake of simplicity we shall disregard the magnetic 
interaction H mag for a moment. Secondly, the s-wave scattering lengths are formally de- 
termined from the so-called T(ransition)-matrix which can be computed perturbatively by 
using the diagrammatic techniques. Moreover, it is known that the energy correction due 
to the two-body interactions can be directly related to the matrix elements of the T-matrix 
0]. Hence, one expects that the desired s-wave scattering lengths can be obtained from the 
calculations of energy corrections. In fact, it is not hard to show that, to the second order, 
our results agree with those obtained by the T-matrix approach. However, as the standard 
second-order perturbation methods are only required in our paper, the calculations can be 
greatly simplified. With these remarks in mind we are motivated to compute the energy 
corrections due to H mt . 

We first introduce a class of 2-particle states defined by 
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K>- 0; <p) = -^>^4,o4,o l vac ) 



(9) 



where (p a are constant parameters and |vac) is the Fock vacuum. Such states can be nor- 
malized by imposing the condition <£* a <£ a = M 2 = 1- Quite clearly, the unperturbed energy 
vanishes in the presence of the state Eq.(|9|). It is easy to show that the first-order energy 
corrections due to H int is 
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Next, we consider the second-order correction for the energy. Now, in view of the explicit 
form of ifi n t given in Eq.(^), the only possible intermediate states are those states of two 
non-condensate particles carrying opposite momenta 



|k, -k; a, (3) = aj^aj k |vac) 



(11) 



with which the unperturbed energy is given by 



(k, -k; a, (3 \H \ k, -k; a, (3) = e k + e_ k = 2e k . 



(12) 



Thus the second-order energy correction due to ifi nt is 



E, 



(2i 1 |(0,0;v9|gi n t|k, -k;q, 
2 k ^ 2e k -0 
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The factor 1/2 in Eq. (|T3|) is inserted in order to avoid the double counting of the momentum 
states. Now, we have 



and hence 
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Obviously, the integral in Eq.(|TJJ) diverges as |k| — >oc. Choosing ip a in such a way that 
\tp 2 \ = 0. yields 



9i = E$V, 



(16) 



indicating that gi is proportional to the first-order energy correction due to the two-particle 
interaction H [nt . At this order, gi is related to the corresponding s-wave scattering length 
a x by 

47rai 
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m 
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Hence, to the second order, at is related to g\ by the following equation 

m v ' 

9 /" d 3 k 1 . . 

= 9 '- 9? /(^f^ (18) 

Here g\ will be referred to as the corrected coupling constant of g%. Writing the original 
coupling gi in terms of the corrected coupling gi, we have at the same order 

n f (J/ k 1 , . 

9l = 9l + 9l 7p^ (19) 

which is equivalent to the results demonstrated in the one-component case ||. Next, we 
consider the corrections of g%. Unlike gi, g 2 can not be isolated directly in the present 
formalism. Instead, we shall consider the sum g\ + g 2 which is nothing but the coupling g n 
associated with the density-density interaction as g% = g n + g s , g 2 = —g s - To this end, we 
may take ip 2 = 1 such that g x + g 2 = E-^V and hence 



= 9i+92- (gl + 2g l92 + 3gl) J ^- (20) 

Subtracting Eq.([0|) from Eq. (|20|) and using Eq.(|l8|) again, we have at this order 

f d 3 k 1 

92 = 9 2 + (2» 1 s 2 + 3sS)/p- f - (21) 
Alternatively, g n and g s are related to the corresponding corrected couplings by 

s "=*» + ^ +2 »«)/(S^ (22) 

/ ~ ~2\ /" d 3 k 1 

It should be noted that the corrected coupling constants, g\ and g 2 , are consistent with the 
one-loop corrections obtained by using the Feynman diagram techniques M. These results 
are actually unaltered in the presence of a constant magnetic field. The point is that the 
two-body interaction term H mt , in fact, commutes with the magnetic term H mag . As a 



consequence, despite that the magnetic interaction would, inevitably, introduce a Zeeman 
energy shift to each hyperfme spin state the total Zeeman energy is conserved in the two- 
particle scattering processes. Based on this point, one can easily check that both Eqs.(|T9D 
and plf ) remain correct. 

We now proceed to calculate the ground state energy with the foregoing results. In the 
standard Bogoliubov approximation the operators a a ^ and replaced by the classical 

number & a VV and <&* a \/V respectively, such that |$| 2 = Nq/V = n represents the density 
of condensate particles. Making these replacements into Eqs.(||) and (||) yields 
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Using Eqs.(p3|) and (|4]) we obtain the effective Hamiltonian 

H c s H con -\- H non 



92 \ -> /»2 t t ^.*2 
+ 7T 2^ V a Ak a /3 -k + $ a «,k«a -k 
^ k^O 
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#non = E ( a l,k^«/3 a /3,k + ~ M^d^ap _ k + T^a^k /? -k) ( 27 ) 
k^O ^ 1 1 J 

are the Hamiltonians for the condensate and non-condensate part, respectively. Here the 
matrix elements are given by 
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£a/3 = 6k<W + iBa^f + 9l |$| 2 + 9l%® a + 2c/ 2 $* $/3 
•Ma/9 = + #2$ 2 <W- 



(28) 



Note that H e s is precisely the Hartree-Fock-Bogoliubov Hamiltonian in the standard Bo- 
goliubov approximation [0] whose ground state structure can be determined by minimizing 
the integrand in Eq. (Pop . As a result, two different ground state structures are found 

(29) 



® =yfm (l/V2,i/V2,0) for n g 2 > 
which is referred to as the "ferromagnetic" state and 

<fr =y/n ~ ( cos 9, i si n 0) for n g 2 < — 
as the "polar" state. Here the cosine and the sine in Eq.(|30"D are given by 



(30) 
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Since if e fr is quadratic in a a k and a a ^ we can diagonalize this Hamiltonian by using the 
generalized Bogoliubov transformation 



where i is the mode index and 



and annihilation operators and 6j c w , the Hamiltonian takes the following form 



V \u {i) b {i) - v {i) b Ki) 



(32) 
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1. In terms of the quasiparticle creation 



i k^O V 



.(*) 



(33) 



We now define the ground state which is annihilated by all b-£ i.e., 6^ |GND) = 0, such that 
the ground state energy is found to be 

d 3 k 



E. 
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where the last integral indicates the energy shift due to the quasiparticle excitations. To 
calculate the ground state energy, one needs to know precisely the values of and v ^ k for 
the quasiparticle modes. This can be done by using the standard Hartree-Fock-Bogoliubov 
mean-field method, and we have calculated E± and for the quasiparticle modes. In the 
following, we devote our attention to the case in which the ground state is "polar", since 
for the "ferromagnetic" case the results are identical the those of the one-component scalar 
BEC || and can be obtained from the "polar" case by setting Q = —11^2. As a result, 
we find that the low lying excitations can be described by two gapless modes £7 k and one 



massive mode E- 



(0). 



E, 



(±> 



E 



(0) 



(35) 



for which the corresponding nonvanishing distribution functions are given by J7| 



where 
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r/ {±) = n - 2g&) + n + 2g 2 ) y/l - (fi/no \g 2 \f 



(40) 



Moreover, with the condensate wavefunctions described in Eqs.([3(J) and ([31]), we obtain the 
following results: 

ft 2 \ „ 



1*1 ="-0, 



I #2 I 



(41) 
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However, one see that all the three integrals 



d 3 k 
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are divergent when |k| — > oo. Note also that they are essentially of second-order in the cou- 



pling constants g\ and g 2 - To eliminate these ultraviolet divergences we substitute Eqs. (|i~9|) 
and (|21|) into Eq.(R4|). The resulting expression for the ground state energy is 
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Note that the condensate wavefunction $ determined by minimizing the sum of terms in the 
first line of Eq.(^) has the same form as that in Eq.([Tj]) except that the corrected coupling 
constants are substituted instead. Furthermore, the last term in Eq.(^2|) can be expanded 
in powers of the corrected coupling constants. Since Eq.(|42"D is valid only up to the second 
order in the corrected couplings, it suffices to substitute gi = gi, g 2 — 92 m the expressions 
for E^ and t>^ k , i.e., in Eqs.(E3)-(H0). On these grounds, we are now ready to calculate 



.Egnd given by Eq. (^2|) . First, we note that 



+ 9 



~g{ + 2~ 9l (g, +g 2 )[l 
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and hence 
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The integral in Eq.(f4"2"D is then equal to 
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Note that the first two terms are the same as that of the one-component case 



d 3 k _ (±)2 / 1 1 



,5/ 2 

m V 15y7r 



128 _ 3 (±) 5 / 2 " 
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where = mg^' /An are the corresponding s-wave scattering wavelengths. The last 
integral can be expressed as 



2 n ° V I 92 ~ 4 



(0) 



Ck + ™0 |<?2 1 
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where £ = fl/no Ifitel* anc ^ the function of integral is defined as 

15^2 
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(48) 



which is a monotonically increasing function that can not be analytically evaluated in gen- 
eral. 

Finally, using Eq.(]£TD we get 
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and therefore the ground state energy is given by 
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2nn 2 V 
m 
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(50) 



where a n = mg n /4:7r, a s = mg s j '4-7T. The terms proportional to a^ 5 ^ 2 are caused by the two 
gapless modes and have the same form of the phonon-like mode in the one-component BEC. 
The last term in Eq.(|50"D is due to the massive mode, which depends solely on the scattering 
length a s for the spin exchange channel and is suppressed by the increasing magnetic field. 
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In conclusion, we have analytically calculated the ground state energy of a homogeneous 
spinor BEC with hyperfine spin / = 1 based on the Bogoliubov approximation. In this weakly 
interacting system, the two-body interactions are described by the hard-core collisions and 
the spin exchange interaction which are characterized by the coupling constants g\ and g 2 
respectively. Using the second-order perturbation methods, the two bare coupling constants 
gi and #2, are expressed in terms of their corrected ones, g\ and g 2 which are directly related 
to the physically measurable s-wave wavelengths for the corresponding scattering channels. 
It is found that the correction of g\ has the same form as that of the one-component scalar 
BEC. However, the correction of g 2 is more complicated and has dependence on the corrected 
coupling constant g x . With the corrected coupling constants, we are able to show that the 
ultraviolet divergence occurring in the calculation of ground state energy can be completely 
removed. 

This work is supported by the National Science Council, Taiwan under Grant NSC-88- 
2112-M-018-004. 
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